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HAMMERSLEY’S HARNESS PROCESS: INVARIANT DISTRIBUTIONS 

AND HEIGHT FLUCTUATIONS 

TIMO SEPPALAINEN AND YUN ZHAI 


Abstract. We study the invariant distributions of Hammersley’s serial harness process 
in all dimensions and height fluctuations in one dimension. Subject to mild moment 
assumptions there is essentially one unique invariant distribution, and all other invariant 
distributions are obtained by adding harmonic functions of the averaging kernel. We 
identify one Gaussian case where the invariant distribution is i.i.d. Height fluctuations in 
one dimension obey the stochastic heat equation with additive noise (Edwards-Wilkinson 
universality). We prove this for correlated initial data subject to polynomial decay of 
strong mixing coefficients, including process-level tightness in the Skorohod space of space- 
time trajectories. 


1. Introduction 

Thinking about the crystalline structure of metals around 1956 led J. M. Hammersley to 
formulate the serial harness. This process ht{x) evolves on the lattice via the equation 

(1-1) ht+i{x) = '^p{y- x)ht{y) + Ct+i{x) 

y 

where p is a symmetric random walk kernel on and are i.i.d. random variables 

with mean zero and hnite variance. These ideas were recorded later in article [15] that 
contains, among other things, calculations that relate the fluctuations of the process to the 
behavior of the random walk kernel. 

Toom [35] studied the convergence of ht{x) as t —)• oo, as a function of the tail of the 
noise ^t(x). [35] also contains several references to the physics literature. 

A natural continuous-time version of the process applies the local update (1.1) at the 
epochs of Poisson processes attached to lattice points x. The ergodic properties of the 
continuous-time process were investigated by Hsiao [17, 18], first for the Gaussian case and 
then more generally. In particular, [18] recorded the order of magnitude of the fluctuations 
of ht{x) from the identically 0 initial condition: in d = 1, -y/log t in d = 2, and bounded 

in d > 3. In d > 3, [18] showed (i) uniqueness of an invariant distribution that is invariant 
under spatial shifts and has given finite mean and variance, and (ii) convergence to such 
equilibrium from initial distributions that are invariant and ergodic under spatial shifts and 
possess a finite second moment. 


Date'. April 28, 2015. 

Key words and phrases. Harness, Gaussian process, Edwards-Wilkinson universality class, random walk, 
fractional Brownian motion, fluctuations, interface, process tightness, strong mixing coefficients, linear pro¬ 
cess, harmonic crystal, stochastic heat equation. 


1 




2 


TIMO SEPPALAINEN AND YUN ZHAI 


In the same continuous-time setting, Ferrari and Niederhauser [13] introduced a dual rep¬ 
resentation of the harness process in terms of backward random walks. Through martingale 
techniques and the random walk representation, [13] proved convergence of the process from 
flat initial profiles both on and on subsets of and obtained some rates of convergence. 
In d S {1,2}, instead of ht itself, [13] considered ht — htifi) (the process as seen from the 
height at the origin) or the pinned process with boundary condition /if(0) = 0. In partic¬ 
ular, [13] identified the invariant distributions of the harness with Gaussian noise in d > 3 
as Gaussian Gibbs fields (harmonic crystals) studied by [6]. 

Our paper works with the discrete time process (1.1), with a general finitely supported 
random walk kernel p. We establish that in one space dimension the height fluctuations 
of the process obey Edwards-Wilkinson universality. This means that on the space and 
time scale n, fluctuations of the height are of order spatial correlations occur on the 
scale and limit distributions are Gaussian. In particular, height fluctuations of the 

stationary process converge to fractional Brownian motion with Hurst parameter 1/4. We 
do not address limits in higher dimensions. 

We prove the height fluctuations for spatially stationary, correlated initial data, subject to 
polynomial decay of strong mixing coefficients. Our proofs use the discrete-time counterpart 
of the dual representation of Ferrari and Niederhauser [13] and employ a GLT for linear 
processes due to Peligrad and Utev [26]. We obtain also process-level tightness in the 
Skorohod space of space-time trajectories. 

As preparation for the fluctuation results we study the invariant distributions of the 
process. In this part there is no advantage in restricting to one dimension so we do it in 
general. The difference with past work is that we consider the increment process, as is 
necessary for invariant distributions in d S {1, 2}. In any case, the increment process is the 
natural object to study because total increment is conserved. (In one dimension, if height 
h{x) jumps, the changes to the left and right increments h[x) — h{x — 1) and h{x -|-1) — h{x) 
cancel each other.) For a given kernel p and noise distribution there is essentially one 
invariant distribution and its transformations by adding harmonic functions. We prove this 
uniqueness under a condition on the growth of the first moment of the increments on the 
lattice Z'^. 

To summarize, these are the contributions of this paper. 

(i) For invariant distributions we cover low dimensions d = {1,2}. Our convergence and 
uniqueness results require weaker moment assumptions than earlier results and we do not 
assume spatial shift invariance for uniqueness. 

(ii) We prove that the harness process obeys EW universality. Eor EW class height 
fluctuations in general, past work is for product form initial distributions, while our paper 
covers some correlated initial data. In particular, we see that fractional Brownian motion 
with Hurst parameter 1/4 arises from a stationary process also when the invariant distri¬ 
bution is not of product form. We also obtain process-level tightness in space and time, 
which was done earlier only for independent walks. 

Past proofs of one-dimensional EW class fluctuations for an interface or a particle cur¬ 
rent covered independent particles [21, 32, 33], independent particles in a static random 
environment [27] and in a dynamic random environment [20], the random average process 
[2, 12], and a recent continuum example from Howitt-Warren flows [37]. In the context of 
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colliding particles where the position of a tagged particle is a function of current, related 
work goes back to Harris [16] and later in [9]. Fluctuations of order and fractional 
Brownian motion limits have also been identified in the simple symmetric exclusion process 
since the classic work of Arratia [1] and later in [8, 19, 25, 30]. 

The EW class should be contrasted with the KPZ (Kardar-Parisi-Zhang) class where the 
exponents are 1/3 for height or current fluctuations and 2/3 for spatial correlations. Recent 
reviews of KPZ appear in [4, 7]. Roughly speaking, EW fluctuations are expected when 
the macroscopic velocity of the interface is a linear function of the slope (as in Theorem 
2.1 below), while KPZ fluctuations are expected when this connection is nonlinear. 

Organization of the paper. Section 2 describes the model and the results on the invari¬ 
ant distributions, and Section 3 contains the height fluctuations in d = 1. The remainder 
of the paper covers the proofs: Section 4 for invariant distributions. Section 5 for finite¬ 
dimensional convergence of height fluctuations, and Section 6 for process-level tightness of 
height fluctuations. 

Notation and conventions. We collect here some items for quick reference. Z_|_ = 
{0,1,2, 3,... }, N = {1,2,3,...}, [n] = {1,2, ...,n}. The d-dimensional integer lattice 
is = {x = (xi,..., Xd) '■ each Xi £ Zj. jxj is the absolute value of x € M or the Euclidean 
norm of x G The imaginary unit is i = ^—1. C is a finite positive constant whose 
value may change from line to line. 

X ~ /i means that random variable X has distribution /r: P{X G A) = /r(A) for Borel 
sets A on the state space of X. A4i is the space of Borel probability measures on 
X is the subspace of measures that are invariant for the increment process defined below by 
equation (2.14). Transition probability p^{x,y) can be written p^ y to save space. 

E represents generic expectation, and and Var^ denote expectation and variance 
under probability measure p. There are three particular expectation operators, (i) E is 
expectation under the probability measure P of the i.i.d. variables ^ = {^t(x)}igg that 
drive the dynamics, (ii) P is the probability measure and E the expectation of the harness 
process. The superscript on and EX identifies the initial increment configuration ry or the 
initial distribution u. P is the ^-marginal of P. (iii) In the proofs of the height fluctuations, 
P with expectation E refers to the random walks XI of the dual representation of the 
harness process. 

2. The harness process and its invariant distributions 

2.1. The model. Fix a dimension d G N. The state of the harness process is a real¬ 
valued height function /i : Z'^ —)■ M that evolves randomly in discrete time. The evolution is 
determined by a fixed weight vector {p(x)} 3 ,ggd and a collection of i.i.d. “noise” variables 
{?i(^)}tez xeZ''- The weight vector {p{x)} is a finitely supported nondegenerate probability 
vector on Z'^ with these properties: 

0 < p{x) < 1, ^ p{x) = 1, 3M < oo such that p{x) = 0 for jxj > M, and 

Vu G Z'^, the smallest additive subgroup of Z'^ that contains the translated 
support {u + X : p{x) > 0} is Z'^ itself. 


( 2 . 1 ) 
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The last property is strong aperiodicity in Spitzer’s terminology [34, p. 42]. In d = 1 it is 
the same as requiring that p have span 1 [10, Sect. 3.5]. 

It is useful to think of p{x) as a random walk transition probability p{x,y) = p{y — x). 
Multistep transitions are denoted by = lx=y, P^ix,y) = p{x,y), and for k>2, 

P^{x,y)= ^ p{x,xi)p{xi,x 2 )---p{xk-i,y). 

Let vi = J2xez,‘^ xp(x) denote the mean (vector) of p. In d = 1 the variance is 

(2.2) af = '^^(x — ui)^p(x), assumed > 0. 

The state of the height process at time f £ Z is denoted by ht = {ht{x)}x£zd' S 
Given ht, the evolution from time t to t + 1 is governed by the equation 

(2.3) ht+i{x)= '^p{x,y)ht{y)+^t+i{x), x £ Z*^. 

yeZ'i 

Since adding a constant to Ct{x) would simply add a constant speed to ht, we assume that 
IE[.^ 4 (x)] = 0. Our results also require (at least) square-integability. We summarize these 
assumptions as 

^ = {Ctix)}t£z x&<^ i-i-d- real-valued mean-zero random variables 

(2.4) „ ’ „ 

with cr| = E[ lCt(x)l ] < oo. 

The probability measure on the variables ^ is denoted by P and expectation by E, while for 
the harness process we write P and E. P is the ^-marginal of P. 

That the harness process should obey EW universality is indicated by a hydrodynamic 
limit described by a linear partial differential equation — ui-Vu = 0. 

Theorem 2.1 (Hydrodynamic limit). Let uq be a continuous function on and define 
u{t,x) = uo{x -I- tvi). Let {hf} be a sequence of harness processes whose initial height 
functions approximate uq locally uniformly in probability: 

lim p| sup |n“^/io([nxj) — uo(x)| > e| = 0 Ve > 0, ii < oo. 

Then 

lim P{ ( \nx \) — u{t, x) | > e| = 0 V e > 0, (f, x) £ M+ x 

We omit the straightforward proof of the theorem (see [38]). 

2.2. Invariant distributions. We focus on the invariant distributions of the increment 
process because those exist in all dimensions. Once formulated, the increment process 
stands on its own and is more general because not every process that obeys the increment 
evolution comes from a height process. We study the invariant distributions of this more 
general “increment process”. At the end of the section we make brief contact with the 
height process. 

Given a height process /i. that obeys (2.3), define increment variables by 

(2.5) r]t{x, y) = hfiy) - hfix) for x, y £ Z'^. 
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Increment variables satisfy the additivity relation 

(2.6) r]t{x,y)+ r]t{y,z) = T]t{x,z) for x, y, z G 
and by (2.3) the Markovian evolution 

(2.7) Vt{x,y) = ^ p( 0 ,z) 7 /j_i(x + 2;,y + z) + Ct{y) - ^t{x). 

Z&lA 


Conversely, if variables r]t{x,y) satisfy additivity (2.6) at f = 0 and (2.7) for f G N, then 
(2.5) holds for the height process h, defined by ho{0) = 0, ho{x) = r]o{0,x), and (2.3). 
Iterate (2.7) backward in time: 

Vt{x,y) = ^ p{0, z)r]t-i{x + z,y + z) + ^t{y) - it{x) 


(2.8) .-1 

= • • • = X] P% 0 ,z)rit-s{x + z,y + z) + -P^{x,z)) 

zelP k=o zezd- 

for X, y G and s G N. Imagine taking s —)• oo with zero initial condition “y_oo = 0” in 
the infinite past. This suggests the definition 


(2.9) At{x,y) = '^'^ ^t-k{z){p’'(.y,z)-p’^{x,z)), x,yeZ^,teZ. 

k=0 zSZ'* 

The series of independent mean zero variables on the right of (2.9) converges almost surely 
and in L^(P) because the sum of variances converges (Theorem 2.5.3 in [10]): 


( 2 . 10 ) 


IE[Ai(x,y)^] = X] (/(y> 2 :) -y^(x,z))^ 

k=0 


= 2c 7| ^(y''(0,0) - y^(x - y, 0)) = 2f7|a(x - y) 
k=0 

where we introduced the symmetric random walk kernel 

(2.11) y(x,y) = y(0,y-x) = ^p( 0 , 2 :)p(x,y + z), x,yeZ'^, 

z&Zd- 


and its potential kernel 

CX) 

(2.12) a(x) = ^[y'^(0,0) - q'^{x, 0)], x G Z'^. 

k=0 

The strong aperiodicity assumption in (2.1) guarantees that Z'^ is the smallest subgroup 
that contains the support {x : q{x) > 0}. Convergence of (2.12) under this property is 
classical: see Section 28 in [34] for d = 1 and Section 12 for d = 2. In d > 3 convergence 
follows from the transience of the g-walk. 
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Remark 2.2. The convergence issue here is real. Suppose p(—1) +p(l) = 1 in d = 1 so 
that q is supported on {—2,0,2} and assumption (2.1) is violated. Then from (2.9) 

OO OO 

(0, 1) = X] Z] 

A:=0 z^Tj 

splits into two independent sums because for a given {k, z) pair at most one of p^(l, z) and 
p^(0,z) is nonzero. That this series cannot converge even in distribution in d = 1 can be 
shown with the characteristic function argument of the proof of Theorem 2.11(a) in Section 

4. 


One checks from the definition (2.9) that the variables At{x,y) satisfy additivity (2.6) 
and the evolution equation 

(2.13) At+i(x, y) = JJp(0, z)At(x + z,y + z) + 6+i(y) - &+i(x), x,y t£ Z. 

Z 

Note that At{x,y) is a function of {^s}s:s<t and thereby independent of ^t+i. Hence (2.13) 
describes a Markovian evolution. {At{x,y)} is an increment process of a height process. 

By additivity (2.6) it is enough to keep track of nearest-neighbor increments. Let us 
simplify the process to yt = (df(^)}a:eZ‘^ with the vector-valued spin variable r]t{x) = 
{rjt{x — Sj, x))jgjrf] € (The notation is [d] = {1,2,..., d}.) The state space of the time-t 
configuration ry is (M'^)^'^. For y — x = Ci the evolution equation (2.7) specializes to 

(2.14) yt+i{x - ei,x) = '^p{0,z)r]t{x + z - ei,x + z) -1- Ct+i{x) - Ct+i{x - e^), 

Z 

for X € 7^^ and i = l,...,d. We now consider the (M'^)^'^-valued process ry defined by 

(2.14) , without imposing the additivity requirement (2.6). 

Let us establish some standard notation and terminology. A generic element of the state 
space (M'^)^ of process (2.14) is denoted hy y = {y{x))x^id with y{x) = (T?(a: — Cj, x))jg[rf] € 
A4i = A4i((M'^)^'^) is the space of probability measures on (M'^)^'^. A measure p G A4i 
is invariant for process (2.14) if equation (2.14) preserves this distribution: that is, liyo ^ p 
and noise = {Ciix)}xezd is independent of yo, then yi defined by (2.14) satisfies again 
r/i ~ /r. I denotes the convex set of probability measures in A4i that are invariant for 
process (2.14). 

Write {0x,t}xezd tez shift mappings in space (x-index) and time (t-index). When only 
space or only time is shifted, the other index is set to zero. For example, for y G (M'^)^'*, 
{6a,oy){x — ei,x) = y{x + a — ei,x + a). A probability measure p G A4i is (spatially) 
shift-invariant if p{6xflA) = p{A) for all Borel sets A C (M'^)^'^ and x G 7^. A Borel set 
B C (M'^)^ is shift-invariant if OxpB = R Vx G Z'^. A shift-invariant probability measure 
p G A4i is ergodic if p{B) G {0,1} for every invariant Borel set B. A shift-invariant 
probability measure p on (M'^)^ is ergodic under the individual space-time shift Oz^s if 
p{A) G {0,1} for every Borel set A that satisfies Oz^sA = A. 

The Markov chain y, defined by (2.14) without the additivity requirement (2.6) gives 
us a larger class of processes than the increment processes coming from height processes. 
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But we find that, under very lenient moment assumptions, all stationary ry. processes are 
obtained by adding harmonic functions to A. of (2.9). 

The hrst theorem summarizes what was developed above. Set 

(2.15) At{x) = {At{x - ei,x))i(,[d] e and A* = {At{x)} G (M'^)^^ 

The invariance and ergodicity claim below follows from At (a:)(0 = Ao(0)(6»a,,f^). Let tto G X 
denote the distribution of A^. 

Theorem 2.3. Assume (2.1) and (2.4). Then the series in (2.9) converges for almost 
every ^ and in L^(P). For almost every ^ the variables At{x,y) satisfy additivity (2.6) 
and the evolution equation (2.13). The process A. = {Af(x)}^^^ defined in (2.15) is 
a stationary version of the Markov chain (2.14) and is invariant and ergodic under each 
individual space-time shift mapping for {x,t) / (0,0). 

Next we add harmonic functions to the variables At{x — ei,x) to produce invariant 
distributions with nonzero means and address uniqueness and convergence. 

A function u : —)• M is harmonic for transition probability kernel p if 

(2.16) p{x,y)v{y) = v{x) for all x G 

Constants are the only bounded harmonic functions (Lemma A.l in the appendix) but there 
can be many unbounded ones. Let TLd denote the space of harmonic functions u ^ 
in other words, the space of vectors u = (tti,... ,Ud) of d real-valued harmonic functions 
Ui ^ M. Given u G TLd, dehne a process Af = u -T At hy 

(2.17) Af{x — ei,x) = Ui{x)-\-At{x — ei,x), t G Z, x G Z'^, i G [d], 

where At(x — ei,x) is defined by (2.9). A“ is another time-stationary and ergodic process 
whose evolution obeys (2.14). Let 7r„ G X denote the distribution of the conhguration A“. 

For nonconstant u, A“ may fail additivity (2.6). We are now considering the broader 
class of processes dehned by (2.14), regardless of whether this process comes from a height 
process. 

We record the invariant distributions vr^ in the next theorem. That tTu is an extreme 
point of X means that if vr^ = 6/i -|- (1 — h)v for 0 < 6 < 1 and /r, p G X, then y, = v = 7r„. 

Theorem 2.4. (Existence.) Assume (2.1) and (2.4). Then the measures {tt^ : u G TLd} ore 
extreme points of I. They are also the unique invariant distributions of minimal variance: 
for any p G X with finite variances, 

(2.18) Var'^[r/(x — ej,x)] > Var^°[77(x — ej,x)], 

and equality holds for all i G [d] and x G Z'^ iff ^ {'^u ■ u G TLd}- 

The next theorem shows that these are the unique extreme invariant distributions under 
a growth condition on the centered first moment of ? 7 (x). 

Theorem 2.5. (Uniqueness.) Assume (2.1) and (2.4). Let p G X satisfy these properties: 
E'^\r]{x)\ < cx) Vx G Z*^ and there exists u G TLd such that 

(2.19) lim • max } \ r]{x) — u{x) \ } = 0. 

r^QO x&7-‘^:\x\<r 
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Then (2.19) holds also if u{x) is replaced by u{x) = E'^[r]{x)]. There exists a probability 
measure 7 on such that v = '^u+a' 7 (doi). 

Above vr^+Q. is the measure vr^ for the harmonic function u(x) = u(x) + a in Tld- 

Remark 2.6. The key to the proof of Theorem 2.5 is the local central limit theorem applied 
to the random walk p*(0, x). The factor in assumption (2.19) comes from a smoothness 
bound of Gamkrelidze [14, Theorem 4]: 

d 

( 2 . 20 ) ^ ® s ^+- 

l=l 

For an irreducible, aperiodic kernel p the same estimate is given in Prop. 2.4.2 of [22]. 

The third item is convergence to an invariant distribution. We currently have a result 
only when the centered initial distribution is of this type: 

C = {C{x)}x£Zd is an (M*^)^ -valued random configuration whose distribution 

(2.21) is invariant and ergodic under the spatial shift group {0x,o}xeZ‘^i 
the M'^-valued variable ({x) = {C{xd))ie[d] h^-s rnean Ef{x) = 0. 

Theorem 2.7. (Convergence.) Assume (2.1) and (2.4). Let {rit}t&z+ be an -valued 

process whose evolution is defined by (2.14). Assume that the initial distribution of the 
process ry is of the following form: r/o(x — e,, x) = ufix) -\-C,{x, i) for x and i = 1,..., d 
where u = {ui,... ,Ud) G Tld and f is as in (2.21). Let pt denote the distribution of rjt- 
Then, as t —)• 00, pt ^ weakly in M.i . 

Remark 2.8 (Spatially invariant case). Let us spell out the most natural special case. 
Given a constant a £ there is a unique tTc G X that is invariant and ergodic for the 
process pt defined by (2.14), invariant and ergodic under the spatial shift group {Ox,o}xezdj 
and has mean E'^°‘ [p{x)] = a. tTq, is the distribution of the (M'^)^ -valued random con¬ 
figuration {C,{x)}x^'Ed defined by f{x) = a (A(x — ej,x))jg[rf]. Furthermore, the process 
(2.14) started with an arbitrary mean-a spatially invariant and ergodic initial distribution 
converges weakly to tTq,. 

In the one-dimensional case that is the subject of the next section, ttq is the distribution 
of the sequence {AQ{x)}xez where Ao(x) = Ao(x — l,x) is defined by the series (2.9). The 
next theorem collects properties of the shift-invariant covariance 

Vo{x,y) = E[Ao(x)Ao(y)] = E^°[p{x)p{y)] = E^°[r]{0)p{y - x)] = Vo{0,y - x) 

of TTo in d = 1. The span 1 assumption on p implies that the support of q cannot lie in 
a proper subgroup of Z, and thereby the characteristic function (pqiO) = 
equals 1 only at 0 G 27rZ. 

Theorem 2.9. Assume (2.1), (2.4), and d = 1. Then 3 constants 0 < A, c < 00 such that 
|Vo(0, x)| < We have the identities 

2 

( 2 . 22 ) = al[a{x-l) + a{x+ 1 )-2a{x)] = — I ^^ x G Z, 
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and 

2 

(2.23) ^Vo(0,x) = ^ 

xez *^1 

where the series above is absolutely convergent. 

The first equality in (2.22) comes from a series like (2.10) and the second equality is from 
p. 355 of Spitzer [34]. The rest is proved in Section 4 from properties of the kernel a{x). 

Remark 2.10 [d = 1 Gaussian case). Presently we cannot say more about the distribution 
TTo, except in the Gaussian case. If {^t(x)} are centered Gaussian variables, then ttq is the 
mean zero Gaussian measure with covariance (2.22). Furthermore, ifp(0, z)+p{t), z + 1) = 1 
for some z ^ "L, then q is supported on { — 1,0,1} and (2.22) shows that Vb(0)2:) = 0 for 
X / 0. In other words, the variables {Ao(x)} are uncorrelated, and thereby in the Gaussian 
case they are i.i.d. 

We turn briefly to the invariant distributions of the height process ht- The key point is 
that in d € {1, 2} there are none. By analogy with (2.9), define 

CX) 

(2.24) Xt{x) = X] X] it-k{z)p''{x,z), tez. 

k=0 

Gomputation of the sum of variances on the right gives 

OO OO 

Hxt{xf] = o-| ^ ^ /(x, zf = q^io, 0) = a|G(0,0) 

k=0 k=0 

where the Green function 

CXD 

(2.25) G{x, y) = ^ q^ix, y), x,y e 

k=0 

converges by transience in d > 3. 

Theorem 2.11. Assume (2.1) and (2.4). Consider the height process h, defined by (2.3). 

(a) In dimensions d = 1 and d = 2 this process has no invariant distributions. 

(b) In dimensions d > 3 the series in (2.24) converges for almost every ^ and in 

L^(P). For almost every ^ the variables Xt{x) satisfy equation (2.3). The process x. = 
{Xtix)}t£z xei'i ® stationary version of the Markov chain (2.3) and is invariant and 
ergodic under each individual space-time shift mapping 9x,t for {x,t) (0,0). 

Part (b) of the theorem is clear. Part (a) is proved in Section 4. 

The stationary height and increment processes (2.24) and (2.9) are obviously connected 
by At(x,y) = xt{y) - Xt{x). In fact any bi-infinite process {'ntix,y)}tei,,x,yei.d- that sat¬ 
isfies additivity (2.6) and evolution (2.7) \/t £ Z comes from a bi-infinite height process 
{ht}tez, uniquely determined once a value /io(0) is chosen: the forward process {ht}tez+ 
is determined by (2.5) and the evolution (2.3). We extend the height process backward in 
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time inductively. Assuming that {hs}s>t has been constructed, extend to time t — 1 by first 
defining the value 

hj_i(0) = ht{z) - '^p{z,y)r]t-i{0,y) - Ct{z) 
y 

where 2 : S is any point, and then ht-i{x) = /if_i(0) + for all x / 0. That the 

definition of /if_i(0) is independent of 2 ; is equivalent to (2.7). That ht now comes from 
ht-i by (2.3) is immediate from the definitions. 

3. Height fluctuations in one dimension 

Restrict to dimension d = 1. Assume that the initial height function ho = {ho{x)}x£Z 
is normalized by ho{0) = 0 and that the distribution of the initial increment configuration 
Vo = {%(2;) = ho{x) — ho{x — l)}xez is invariant under translations of the spatial index x. 
The mean, variance, and series of covariances of the initial increment variables are 

(3.1) /To = E[r/o(x)], o-Q = Var[? 7 o(x)], and ^ Cov[?7o(0), r/o(x)]. 

Our assumptions will guarantee absolute convergence of the series. Then go > 0 because it 
is the limit of n“^Var[r7o(l) + • • • + ?/o(^)]- 

Let b = —vi = — xp(0, x). Scale space by ^/n and time by n and consider the scaled 

and centered space-time height process 

(3.2) hn{t,r)=n~^/'^ {h\nt\{[ry/n\ + [ntb\) - yory/n] , {t,r)eR+xR. 

We prove that this process converges to a Gaussian process. 

The limit process has three natural descriptions: as a Gaussian process, as a sum of two 
stochastic integrals, and as the solution of the stochastic heat equation with additive noise. 
The proof is based on the Gaussian process description. Denote the centered Gaussian p.d.f 
and c.d.f with variance by 


(3.3) 


nx 

= / v’u^{y)dy, 

J — 00 

and let 

(3.4) 

Tj ,2 (x) = (x) - X (1 - ^^2 (x)), 

X G M. 


Define two positive definite functions Ti and r 2 for {s,q), {t,r) G M+ x M by 
(3-5) ri((s,g), (t,r)) = T^ 2 (i+^)(r - q) - T^ 2 |i_^|(r - q) 

and 

(3-6) r2((s,g),(t,r)) = T^2^(-g)-h T^2i(r) - T^2(i+^)(r - g). 

Let {Z{t,r) : {t,r) G M+ x M} be the mean zero Gaussian process with covariance 

2 

(3.7) E[Z{s,q)Z{t,r)] = ^Ti{{t,r),{s,q)) + goT 2 {{t,r),{s,q)). 
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Equivalently, this process is given by 

(3.8) Z{t,r) = ^ jj - x)dW{s,x) + J ip„2^{r 

[0,t]xR K 


x)B{x)dx, 


where {W{t,r) : t G M_|_,r G M} is a two-parameter Brownian motion and {B{r) : r G M} a 
two-sided Brownian motion, and W and B are independent. In the fluctuation limit of 
described by (3.7) and (3.8), the first term comes from the fluctuations of ^ (the dynamics), 
and the second term comes from the fluctuations of rjo propagated by the dynamics. 

Z is also the unique mild solution [36] of the stochastic partial differential equation 

(3.9) onM+xM, Z(0,r) = y/^B(r). 

We consider three different hypotheses on the initial increments {? 7 o(x)}: (a) i.i.d., (b) 
strongly mixing, and (c) the invariant distribution tto of Theorem 2.3, defined by 


OO 

(3.10) %(a^) = xeZ. 

yEZ k=0 


Throughout the text these are referred to as cases (a), (b) and (c). In case (c) parameter 
qq disappears from the limit because ^>o = ^\l^\ (2.23). 

Given two sub-u-algebras A and B on a probability space let 


(3.11) a{A,l3)= sup \P{AnB) - P{A)P{B)\. 

A&A,B&B 


For the initial increment sequence {r]{x)} define cr-algebras = o'{r]{x) : m < x < n}, 
and then the strong mixing coefficients 


(3.12) 


ain) =supa{P^_^,^, n G Z+. 

k 


The sequence {r]{x)} is strongly mixing if a{n) —)• 0 as n —)• oo. See Bradley [5] for properties 
of these and other mixing coefficients. 

Here is the convergence of finite-dimensional distributions. 


Theorem 3.1. Assume d = 1, (2.1), (2.4), and E[^t(x)^] < oo. Assume that the initial 
increment sequence {'i]o(x)}xez satisfies one of the assumptions (a), (b), or (c): 

(a) {'rio{x)}x&z is an i.i.d. sequence and E[r/o(ic)^] < oo. 

(b) {'i]oix)}x&z is a strongly mixing stationary sequence and 3(5 > 0 such that E[ |??o(0)P''''^ ] < 
oo and the strong mixing coefficients in (3.12) satisfy 

OO 

(3.13) ^(j -h l)^/^a(j) < oo. 

j=0 

(c) {voixjjxez has the distribution tto of the sequence in (3.10). 

Then, go = Ylxez^^'^idoW^'noix)] > 0 is absolutely convergent. As n ^ oo the finite¬ 
dimensional distributions of the process h^ of (3.2) converge weakly to those of the mean- 
zero Gaussian process Z with covariance (3.7). 
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The convergence in the theorem means that for any fixed G N and (t 2 ,r 2 ), ■ ■ ■, 

(tN,rN) € X 1^) the weak convergence of M'^-valued vectors holds: as n —>■ oo, 

(3.14) {hniti,ri),hn{t2,r2), • • .,hn{tN,rN)) {Z{ti,ri),Z{t2,r2),.. .,Z{tN,rN)) ■ 

Remark 3.2. Case (a) is stated only for completeness, we do not prove it. A natural 
question is whether case (c) is actually covered by case (b). The answer is affirmative in 
the Gaussian case where, Vn G N, a{k) = 0{\k\~"') as |A:| —)• oo. See [38] for details. 

Remark 3.3 (Fractional Brownian motion). In case (c) of the theorem, or in general when¬ 
ever go = the limit of the time-indexed process hn{t, 0) is the Gaussian process with 

covariance 

2 

(3.15) E[Z{s, 0)Z{t, 0)] = + t'/' -\t- ) 

\/27rcrf 

which is a fractional Brownian motion with Hurst parameter 1/4. 

By strengthening the assumptions we upgrade the weak convergence hn ^ Z to process 
level on a compact time-space rectangle Q = [0,T] x [—R, Rj. Z is a continuous process on 
Q. The paths of hn lie in the space D 2 of 2-parameter cadlag paths, defined precisely as 
follows. Given (toi^o) £ Q, define four quadrants by 

Q(to,ro) = {(^> r) e Q :t>to,r > ro}, Qlto,ro) = {(^> r) e Q : t > to,r < ro}, 

Qfto,ro) = {(^> r) € Q :t <to,r < ro}, Qfto,ro) = {(^> r) € Q : t < to,r > ro}. 

Then the path space is dehned by 

D 2 = I f : Q ^ '■ V(to) ’"o) G Q, li™ /(G ^) exists for i G {1, 2, 3,4} 

and liin f{t,r) = f{to,ro)\. 

(‘-'■)6Q(to,ro) ^ 

D 2 is separable and topologically complete under a Skorohod-type metric 

d{f,g) = infmaxdl/ - 50 A||oo , ||A - id||oo), f,g € D 2 , 

where A(t,r) = (Ai(t), A 2 (r)) for strictly increasing, continuous bijections Ai and A 2 . For 
details we refer to Bickel and Wichura [3]. 

Theorem 3.4. Assume d = 1, (2.1), (2.4), and E[^i(a;)^^] < 00 . Assume that the initial 
increment sequence {goix)}xei, satisfies one of the assumptions (a), (b), or (c): 

(a) {go{x)}x& is an i.i.d. sequence and E[t/o(x)^^] < 00 . 

(b) {go{x)}xez is a strongly mixing stationary sequence, and 35 > 0 such that E[ |7?o(0)d^^'^ 
00 and the strong mixing coefficients in (3.12) satisfy 

OO 

(3.16) + l)^°'^^^^/‘^a(j) < 00 . 

j=0 

(c) {go{x)}xe'Z das the distribution tto of the sequence in (3.10). 
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Then on any rectangle Q = [0, T] x [—ii, ii], hn => Z on path space D 2 . 

This completes the description of results and we turn to proofs. 

4. Proofs for invariant distributions 
For X e i € [d], and s < t in Z let 

t 

- ei,x) = ^ ^ik{y)\p^~^{x,y) - p^~'^{x - ei,y)]. 
fc=s+l y 

Then A_oo,t is what we denoted by At in (2.9). As t — s ^ 00 , the random process above 
indexed by {x,i) converges weakly to the configuration Aq ~ vro, by the convergence in 
(2.9). The basic evolution equation (2.14) can be written as 

(4.1) Vt{x - ei,x) = '^p\x,y)r]o{y - ei,y) + Ao,t{x - ei,x), t>0. 

y 

The proof of Theorem 2.4 is contained in the next lemma. 

Lemma 4.1. We have these properties of invariant distributions. 

(a) Let u £ I and assume that under v each variable r]{x — ei,x) has finite mean and 
variance. Then 

(4.2) Var^[r/(x — Cj, x)] > Var’^°[r?(x — Cj, x)]. 

Equality in (4.2) holds for all i € [d] and x iff n £ {tTu ■ u £ TLd]. 

(b) For each u £ Tid, 'Xu is an extreme point ofX. 

Proof, (a) When the process r/t ~ u in (4.1) is stationary, taking expectations on both sides 
shows that Ui{x) = E''[r]{x — ei,x)] is harmonic. Then also the process r/t(x — ei,x) = 
rjfix — ei,x) — Ui{x) satishes (4.1) and is stationary. By stationarity, by the independence 
of the terms on the right of (4.1), by time-shift-invariance of by the L^(P)-convergence 
in (2.9), and finally by the definition of vro: 

\aW[r]{x - ei,x)] = B'"[fit{x - ei,x)^] 

2- 


(4.3) 


= E"' 


'^p\x,y)rifiy - ei,y) 


+ E[Ao,t(x - ei,x)^ 


> E[A_t_o(3: - ei,x)^] 


t —>-00 


E[Ao(x - ej,x)^] = Var""°[?7(x - ei,x)]. 


This gives the inequality claimed in part (a). 

Pick xi ,...,xat £ Tfi, ii,...,iN £ [d], and ai ,...,£ M, and repeat (4.1) for linear 
combinations: 

N N 

'^aiftixi - ei^,Xi) = '^ai'^p\xe,y)fio{y - ei^,y) 

1=1 l=l y 

N 

+ '^aiAodixe-ei^,xe) = St +Df. 

£=l 


(4.4) 
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Equality in (4.2) Vi £ [d],x € and (4.3) imply that —)■ 0. Consequently, on the 

right-hand side of (4.4), the first sum vanishes and the second sum converges in distribution 
to a£Ao(xf — Xf). This implies that under v, the configuration f] has the distribution 
of the configuration Aq, which says that u is among the distributions {vr^ : u G 'H^}- 

(b) To get a contradiction, suppose vr^ = + (1 — /3)i^° for /3 G (0,1) and G I. 

Let rj^ denote the configuration under Let J be a Bernoulli(/3) variable: P{J = 1) = 
/3 = 1 — P{J = 0). Then u = + (1 — /3)p^ is the distribution of the configuration rj'^, 

and by the definition of vr^, the assumption 7r„ = p implies the distributional equality 

(4.5) Ui{x) + Ao{x — ei,x) = rj'^{x — ei,x), x G z G [d]. 

It follows that 1 ]^ and have finite means uf{x) = E’^^[rj^{x — ei,x)] and variances, and 

(4.6) Ui{x) = E''\q\x - ej,x)] = I3u]{x) + (1 - /3)n°(x). 

From (4.5), by expanding the variance, and then by (4.2), 

Var[Ao(x - ej,x)] = £''^[(r/'^(x - ej,x) - Ui{x))^] 

= I3E'^^ [{if-{x - ei,x) - uj{x) + (1 - li){u]{x) - u°(x)})^] 

+ (1 -/3)E‘^°[(v°(a;-ei,x) -u°(x) -/3{u,^(x) -u°(x)})^] 

= /3Var''^[r/^(x - ei,x)] + (1 - /3) Var''° [? 7 °(x - ei,x)] -|-/3(1 - /3){u-(x) - u°(x)}^ 

> Var[Ao(x - ei,x)] +/3(1 - ^){u\{x) - ^“(x)}^. 

This forces hrst = u^, so by (4.6) u = = u^. Second, we get equality of the variances 

Var*^ [r]^{x — ei,x)] = Var*^ [v^ix — ej,x)] = Var[Ao(x — ei,x)] 
which by part (a) forces which equals tt^. 

□ 


We prove Theorem 2.7 next because it will be helpful for a later proof. 

Proof of Theorem 2.7. By (4.1), the assumption on rjo, and harmonicity, 

r]t{x - ei,x) = '^p\x,y)ui{y) '^p\x,y)C{y,i) + Ao,t(x-ei,x) 

(4.7) " 

= Ui{x) + (x,2/)C(y,^) + Ao,i(x - ei,x). 

y 

The terms on the right are independent. As f —>■ 00 , the second term on the right converges 
in Lf to the constant E^(0, z) = 0 by Lemma A.2 from the appendix. As a process indexed 
by (x,z), the last term converges in distribution to Aq. Since the configuration {zzi(x) + 
Ao(x — Cj, x)}jg[^ has distribution vr^, we have shown that pt converges weakly to 
7r„. □ 

For a probability measure p G A4i, let 

A(z/,r)= max F^|? 7 (x)|. 

|a:|<r 
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For configurations r] G let 

(4.8) gt{v,x,i) = '^p\x,y)r]{y - ei,y). 

y 

Lemma 4.2. Let v G A4i satisfy r~^/‘^A{i',r) —)• 0 as r —>■ oo. Then for x G and 
hj e [d], 

lim E''\gt{r],x,i) - gt{v,x -ej,i)\= 0. 

t^oo 

Proof. By Theorem 4 of Gamkrelidze [14] (see also the next to last paragraph of [14]), under 
assumptions (2.1) on the kernel, there is a constant C < oo such that 

d 

(4.9) ^ ® yt G Z+. 

By (2.1) kernel p has finite range M. So by (4.9), 

E^\gt{g,x,i) - gt{g,x - ej,i)\ < ^ \p\0,y) - p\0,y + ej)\ ■ E’^\r]{x + y)\ 

y:\y\<Mt+l 

< Ct~^/‘^A{u, Mt + \x\ + 1) 

and by the assumption the last quantity vanishes as t —)• oo. □ 

For a G let 7r„ = tTu for the constant function u{x) = a. 

Lemma 4.3. Let i/ G X satisfy r~^^‘^A{i',r) —)• 0 as r ^ oo. Then there exists a probability 
measure 7 on such that p = / TTalido). 

Proof. Using (4.8) write the stationary evolution (4.1) for r/t ~ p as 

(4.10) gtix - ei,x) = gt{po, x, i) + Ao,t(x - e*, x). 

Pick lattice points Xk G indices ik G [d], and reals ak, and consider the distributions 
of the linear combinations 

N N N 

^ ^ (^kPtiXk Cij,, Xk) ^ ^ (l/O) Xkj ik) T ^ ^ Otk^0,t{xk ) Xk) S-t Df 

k=l k=l k=l 

and the same thing shifted by —ej: 

N N 

^ ^ (^kPtjXk Cj (^i]^iXk Cj) ^ ^jAk) 

k=l k=l 

N _ _ 

^ ^ CK/j; Ao,t(^/c ; ^k — Sf Df. 

k=l 
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Compare v and its shift by —e,- through characteristic functions. Use time invariance, 


inequality, and ]e''“ 

< 

1 and — e''^ 

lA 

W 

1 

b\ for real 

a, b. 


'[e'-EfcLi akvi^k- 


^'=)] - E'^[e"^^ 

=1 ctkVi^k 


■eP] 1 

= 1 


Ufc- 

_ E'^[e' 

T,k=l OkVt{xk-ej-ei^ 

,^k-ej)j 1 

= 1 

E‘'[e‘'^‘]E[e''^‘ 

]- 

E^[e'''^*]E[e‘-^*] 

1 



< 1 

E^[e‘'^‘] 


] 1 + 1 E[e‘^*] - 

E[e'-'°*] 




< - St 


I E[e'"°‘] - E[e' 




The last line above vanishes as t —)• oo by Lemma 4.2 and because process Ao,t converges 
weakly to process Aq which is invariant under spatial shifts. We have now shown v invariant 
under spatial shifts. 

Let = f fJ,r(d/ 2 ) be the ergodic decomposition of p. L is a probability measure sup¬ 
ported on probability measures // G that are invariant and ergodic under the spatial 
shift group and that have a finite mean a(/i) = E^[rj{x)\. 

Let / be a bounded continuous function on By invariance of v and Theorem 2.7, 

I fdu = E’^[f{rjt)] = I Eq/(r?t)]r(d/r) J (/) r(d/i) = J E^-{f)^{da) 

where 7 is the distribution of the mean of ^ under V{d^). □ 

Proof of Theorem 2.5. By harmonicity of n, the process rjt{x — Cj, x) = rjtix — ei,x) — Ui{x) 
is also time-stationary. By assumption (2.19) its marginal distribution z? G X satisfies the 
growth bound r~^^‘^A{u, r) —?■ 0. Lemma 4.3 applied to v gives the distributional identity 

77(x) = X + Ao(x) where rj ^ u and A G is a 7 -distributed random vector independent 
of Aq. Consequently u{x) = E'^[r]{x)] = u{x) -|- E'^[X\, and the claim about substituting 
u{x) into (2.19) follows. 

An application of Lemma 4.3 to the distribution v of the configuration fit{x) = r]t{x)—u{x) 

gives a random T G such that rj ^ u satisfies {ri{x)}xez<^ = {u{x) + Y + Ao(x)}^^^d. 
With y ~ 7 , this is the same as p = J" nti+alidce). □ 

Proof of Theorem 2.9. Identity (2.22) says that Vo(0)3^) is a Fourier coefficient of (a con¬ 
stant multiple of) the function f{6) = Extend / from [— 7 r, 7 r] to a 

meromorphic function f{z) = complex plane. The only pole in some neigh¬ 
borhood of [—vr, tt] is z = 0. Expansions give f{z) = and we see that z = 0 is 

(7-^Z ) 

a removable singularity. Consequently / is analytic on [—vr, vr] and its Fourier coefficients 
decay exponentially (Prop. 1.2.20 on p. 20 of [28]). 

Identity (2.23) now follows from 


(4.11) 


k 

lim \a{x + k) — a{x)] =±—^ 

2cJ^ 


(P29.2 on p. 354 of [34]). Note that the variance of the g-kernel is 2af (2.2). 


□ 
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Proof of part (a) of Theorem 2.11. By square-integrability (2.4) and Theorem 3.3.8 in [10], 
the characteristic function y?(a) = satisfies 

(p{a) = 1 — + o(a^) as a —>■ 0. 

Suppose {ht}t£Z+ is a time-stationary M^‘*-valued Markov chain that satisfies (2.3). It¬ 
erating (2.3) gives 


(4.12) 


t-i 

ht{x) = 

j/eZ'* fc=o 


On the right-hand side the initial profile /iq is independent of the i.i.d. variables {f,j{y)}. By 
this independence, by the fact that characteristic functions are bounded in absolute value 
by 1, and by the invariance ho = ht, we have, for any 0 < s < t, 

t-i 

|E[e^“^oO)] I < | = \(p{ap^ix,y))\ 

k=s j/gZ'* 

k=S y£jd. 


t-l S. ^ t-1 

-iaV|(l -ho(l)) ^ ^ {p^{x,y)Y\ < expi -ca^^g^(0,0) 

k=s jygZ'^ ^ ^ k=S 

with a constant c > 0. The second equality above is justified by fixing s large enough 
and by the dimension-independent uniform bound p^{x,y) < for k > s (P7.6 on 

p. 72 in [34]). In dimensions d G {1,2} the g-walk is recurrent (T8.1 on p. 83 in [34]), and 
consequently taking t —)> oo above gives E[e''"^°^*^] = 0 for a ^ 0. This contradiction with 
the continuity of a characteristic function at a = 0 shows that in d G {1,2} there can exist 
no time-stationary height process. □ 



5. Height fluctuations: limits of finite-dimensional distributions 

Let {Xl}t£z+ denote a random walk on Z with initial point Xq = i and transition 
probability p{x,y) of (2.1), and let {l}*}tgz+ similarly denote a random walk on Z that 
uses transition q{x,y) of (2.11). Equivalently, Yf = XI — X^ for two independent p-walks 
XI and X^. Under assumption (2.1), the g'-walk also has span 1. 

Probabilities and expectations of these walks are denoted by P and E, always taken 
under fixed rjo and In particular, we have the following “dual” representation of the 
harness process: for t G Z+ and i G Z, 

t 

( 5 . 1 ) ht{i) = E[ho{Xl)]+Y,E[UXU)\: 

k=l 

where we emphasize that ho and ^ have not been averaged over on the right because the E 
acts only on the walk X]. 
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This lemma is a consequence of the local CLT (Theorem 3.5.2 in [10]) and will be used 
several times in the sequel. 


Lemma 5.1. For a mean 0, span 1 random walk Sn on Z with finite variance a £ M, 
and points £ Z such that lim^^oo On /have 


lim —= ^ P{Sk = On) = — / 
1^00 ^ (T^ Jo 


(5.2) 


V2 


: exp 


TTi; 


-—)dv. 

2v) 


The analysis of height fluctuations begins with a decomposition of the scaled height 
function as 


(5.3) hn{t, r) = noHn{t, r) + Fn{t, r) + Sn{t, r) 
where, with y(n) = [ntb\ + [r^/n\, 

(5.4) Hn{t,r) = - ry^), 

[nt\ [ntj 

(5.5) F„{t, r) = n-V- = x) = n"'/- , 

k=l xGIa k=l 

(5.6) Sn{t,r) = j;(r?o(i) - fio) > i) - < *)} • 

iez 

Recall that b = —vi = —E[Xi]. (5.3) follows from the random walk representation 


(5.7) 


hn{t,r)=n 


[nt\ 


k=l 



from (5.1), ho{0) = 0, r]o{k) = ho{k) — ho{k — 1), and from 


- fJ^orVn 


= E 


vvi”) 


1 r -Lxyin) 


> 0 } 


E 


»o(i) l,x»«<o) 


Y1 


- yorVn 


= > i) - ^?7o(i)P(X‘'^”^ <i)- yorVn 

i>0 i<0 


= n^PyoHn{t, r) + n^/^5n(t, r). 


The three terms in (5.3) will be treated separately both for the convergence of finite¬ 
dimensional distributions (this section) and process-level tightness (Section 6). 

Begin by observing that the first term yoF[n{t, r) on the right of (5.3) is irrelevant: since 
E{X^^j) = + y{n) = —6[ntJ -|- [ntb\ -|- [r-^nj, we have Hn{t,r) = 0(n~^P) and 

lim^^oo Aio-f^n(L ^) = 0) uniformly over (t,r). 

Next note that Fn and Sn are independent. Fn depends only on ^ and, in the limit, 
furnishes the first term in (3.7) and (3.8). 5„ depends only on r]o and gives in the limit 
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the second term in (3.7) and (3.8). The weak limits of and are treated separately in 
Propositions 5.2 and 5.3. 


We begin with F^- Let {F{t,r) : t G K_|_,r G M} be a mean-zero Ganssian process with 
covariance 


(5.8) 


E[F{t,r)F{s,q)] = -^Ti{{t,r),(s,q)) 


Here are alternative expressions for Pi [33, Chapter 2]: 


(5.9) 


/ CXD 

<q- x)P{B^ 2 ^ >r-x) 

-oo 

- P{.B 2 ^ <q-x, B 2 ^ >r-x)] dx, 


where Bt is a standard 1-dimensional Brownian motion, and 


(5-10) Ti{{s,q),{t,r)) = ^ J 


1 /-o-ih+s) I I 


: ex' 


ul\t-s\ y/2'Kv ^ '■ 2v 




Proposition 5.2. Assume (2.4) and E[^t(x)^] < oo. Then, as n ^ oo, the finite¬ 
dimensional distributions of the process Fn converge weakly to those of F. 

Proof. First we argue the convergence of the covariance of F^. Let 

^lnsb\-\-lq^\ independent random walks with transition probability p (even 

if (t, r) and {s,q) should happen to coincide). 

For s = t and r, <7 G M, set Xn = [r^/n\ — [qy/n\. Then, 


[ntj 

k=l xGZ 

[nt\ [nt\ — l 

= n-'/V|^(/LntJ-fc(^^^0)=n-VV| ^ qH0,Xn) 

k=l k=0 

^2 


a. 


2 f2afi 1 


n^QO 2(T 


1 Jo 




: exp 


vrn 


{r - gfi 

2v 


\dv = -^Ti{{t,q),{t,r)). 
J CJ^ 


The limit came from (5.2). 
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Next take s < t and set ~ “ IqV^\ ■ By the Markov property, 


[nsj 


E[F^{t,r)Fnis,q)] = ^ P{X^ 


lntb\ + lr^n} _ j^\nsb\ + [qFn\ _ 


-k 


[nsj —k 


= 0 ) 


k=l 


= n-^/^alE 


r L^sJ 


E pt v\ntb\ + \rFn\ _ j^\nsb\ + 'iqFn\ _ „ | ^\ntb\ + \rFn\\ 
\^[nt\—k [ns\—k I [rxtj —[nsj / 




k=l 
r L^sJ 


■- k=l 
[nsj —1 

= n-VV| ^ E[q^{0,X'J]. 

k=0 




We use again Lemma 5.1 to derive the limit. By the CLT, n ^1‘^X'^ ^ + {r — q) 


By a standard construction (Theorem 3.2.2 in [10]), we can find random variables X^ = X'^ 
such that n~^l‘^Xn —)• + {r — q) almost surely. Then by (5.2), 


[nsj —1 


^ 1 r2afs 1 

lim gfc(o,x„)= / ^^expj 

2cjf Jo y/^TTV 1 


(B^2|i_,| + r-g)2 
2n 


} 


dv a.s. 


By the uniform bound q^{0,y) < Ck we can apply dominated convergence. 


lim E[Fnit,r)Fn{s,q)] = lim n ^ q^{0,Xn) 


[nsj — 1 


■- k=0 


(5.11) 


= alE 

2a 




2<7\s 


: exp 


{- 


(B^2|i,| + r-g)2 


2 ralit+s) 

2|t_s| y/27rv 


1 A I 


y/2'Kv ^ 2n 

(r - g)2 


} 


dv 


^ 2v ^ ^ri((t,r), (s,g)) 


The second last equality comes from a convolution of Gaussians. 

Having taken care of the covariance, we turn to derive the weak limit. Fix iV G N, 
(tjjTj) G M+ X M and 9j G ffi for j = 1,..., A^, and arrange the indices so that 0 = to E 
ti < t 2 < ■ ■ ■ < t]sf. Then, with i(k) = t iff [ntjj + 1 < A: < [ntj+ij. 


N N [ntji 

j=l j=l k=l xGX 

N-l L’^i^+iJ N 

= Y1 {[ntjb\+ [rjy/Jl\,x) 

i=0 k=lnt£j-\-l fcGZ 


n 
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[ntjvJ N 

= ^ {[ntjb] + lrjy/n\,x) = ^ T4,fc 

k=l j=l[k)+l xGZ k=l 

where the last equality dehnes variables Vn,k, independent for a fixed re. We apply the 
Lindeberg-Feller theorem (Theorem 3.4.5 in [10]) to {Vn,k}- (5.11) gives 

NatJ N ^2 

(5.12) hm ^ EVl, = Y e^e,4T^{{U,n),itJ,r,)). 

n.—^ ^ ^ /T^ ' 


fc=i 


*.i=i "^1 


As preparation for the second condition of Lindeberg-Feller (the negligibility condition 
(ii) of Theorem 3.4.5 on p. 129 in [10]) we estimate a moment. Abbreviate Pj{x) = 
^\ntj\-k -)- Yrjy/n\^x). Use E[^„(A:)] = 0, E[^t(x)‘^] < oo and q^{x,y) < . 

C ^ \/ 

n\Vn,k\^] < - E ® 

j=£[k)+l xSZ 
N 


(5.13) 


C 


- n Y1 ( E + 


j=£{k)+l x,yez 
N 


xGX 


N 


< ^ E F”‘^J-ho.o)]^ < ^ E 

j=£(fc)+l j=£(fc)-hW‘- ^ 


By Holder’s and Chebyshev’s inequalities and by (5.13), for e > 0, 

C 1 

E[|Fn,fcrj < ■ 

k=l 

0 as re —>■ oo. 


[nlivj [ntff\ 

Y^KkH\Vn,k\>^}] < E 

(5.14) k=i 


N [ntji 

EE 


^ -k)vi 


< 1^(1 +logn) 


ree^ 


(5.12) and (5.14) verify the hypotheses of the Lindeberg-Feller theorem, which implies 
N [nijvJ N 

Tj) = Y^ ^n,k ^ E! 

j=l k=l j=l 

The proof for Proposition 5.2 is complete. □ 

We turn to Sn- The hypotheses on the initial increments are now relevant. Let {S{t,r) : 
t G M+,r G M} be the mean-zero Gaussian process with covariance 

(5.15) F;[5(t,r)5(s,g)] = QoT 2 {it,r), {s, q)). 

The function F 2 has this alternative formulation [33, Chapter 2]: 

/■O 

r 2 ((s, q), {t, r)) = / P{B^ 2 ^ > q- x)P{B„ 2 ^ > r - x) dx 


poo 

+ <Q- x)PiB^ 2 ^ <r-x)dx 


(5.16) 
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where Bt is a standard 1-dimensional Brownian motion. 


Proposition 5.3. Under the conditions in Theorem 3.1, as n ^ oo, the finite-dimensional 
distributions of the process Sn converge weakly to those of S. 

Proof. We have three cases in Theorem 3.1. 

Case (a). The proof of the lemma under the i.i.d. and second moment assumption goes 
via the Lindeberg-Feller theorem and can be found in [38]. 

Case (b). Now assume that the initial increments rjo = {rio{x)}xez are a stationary, 
strongly mixing sequence such that, for some <5 > 0, E|?7o(0)P'’''^ < oo and the strong 
mixing coefficients {a:(j)} of r]Q satisfy < oo. 

We show the distributional convergence 


(5.17) 


N N 

OjSn (tj ,rj)^Y ^ ‘ 

i=i j=i 


for a fixed vector of time-space points {{tj,rj)}i<j<N S (I^+ x and a fixed real vector 
6 = {9i ,..., 9]\f). Rewrite the liner combination as 


N 


(5.18) 
where 

(5.19) 


Y^i^n{tj,rj) = - ho) 

1=1 iez 


^n,i 


= n-V4 1 


N 

^{i>0}Y^jP{X 

1=1 


[ntj bj -I- \rj y/n\ 
\ntj\ 


N 


>i) 

[ntjb} + [rjUn} 




<i) 


1=1 


Our first task is the n —)• oo limit of the variance 

Var[^a„ = Y ^n,jan,kCov{r]o{j),r]o{k)) 


— 

(j„ — 


(5.20) 


YCovipoiO),po{i))Y (^n,kOin,l+k • 




kez 


The next lemma is part of Theorem 1.1 in Rio [29]. The quantile function Qx of |X| is 
defined by 

Qx{u) = inf{a; G M+ : P(|-T| > x) < u}, 0 < u < 1. 

Lemma 5.4. LetX, Y, and XY he integrable random variables and let a = a{a{X),a{Y)). 
Then 


pCX. 

Cov(X,y)l <4 / Qx{u)QY{v)du. 

Jo 


(5.21) 
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We apply this lemma to show that the series of covariances Cov[ryo(0), ??o(^)] is 

absolutely convergent. Let Qniu) denote the quantile function of |r/o(0) — fio\. Then 

ra(\l\) rl 

^|Cov(r?o(0),r?o(^))| / Qrj{u)‘^du < 4 L E ^{u<a{\£\)}Qr]i'^) 




< 4 


/‘(S 

eel 


L{«<a(K|)} 


{ 2 +S)/S 


du 


1 V(2+<5) ^ 1 

/ Qr}{u) 
Jo 


-I 2/(2+5) 


dn 


Since a(n) \ 0 as n —)• oo, we have 

X] l{«<a(hl)}) du = ^ (2j + 1)(2+ )/ dt 

^/7^r77 ^ + 


(5.22) 


OO 


j=0 


By summation by parts, by a{n) > 0, and by the summability assumption on {a(j)}, 

n 


j=0 

n 

= a(0) - (2n + + 1) + X] - {2j - l)(2+<5)/5 

n oo 

<a(0)+C^(i + < a(0) + C ^{j + < oo. 

i=i j=i 

Consequently the quantity in (5.22) is finite. 

For a uniform variable U on (0,1), Qn{U) = |%(0) — fJ-o\ and so 

[ {Qr,{u)f'^^ du = E[\r]o{0) - < oo. 

Jo 

We have shown that 

X]|Cov(r?o(0),r?o(^))| 

£ez 

(5.23) . oo ^ s/{ 2 +s) 

^ i=o ^ 

Next a lemma for the other series in (5.20). 

Lemma 5.5. For a finite constant C 


E[|r?o(0)-//o|'+'']) 


2/(2+5) 


< OO. 


(5.24) 


sup 

neN,fcez 


^ ^ (^n,ie^n,i+k 
ieZ 


< C < oo. 
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For all /c G Z, 

(5.25) lim J2(^n,ian,i+k = ih2^rj2)). 

iez i<ji,j2<N 

Proof. By the Schwarz inequality 

®n,iOn,i+fc ^ Q.n,i 

ieZ i&Z 

and so (5.24) follows from the finite limit in (5.25) for A: = 0. 

To prove (5.25) expand the sum: 

(5.26) 

^ ^ 0‘n,iOtn,i+k 
ieZ 


lGilii2<A' i>0,i+fc>0 


E p,^[ntj^b} + [rj^^} ^ Ap/^Ki2^J + Tj2V^J 


> i + /c) 




E 0n»n E 


IGiU^GA" i<0,i+fc>0 

-1/2 a. a. sr^ + 




+ n 


E 






The middle two terms above are 0{n ^/^) and hence vanish as n —)• oo. 

The individual probabilities converge by the central limit theorem; with i = [xy/n \, 


P{X 


[ntb] + [p\/nj 
[nt\ 


>i) = p 


f^lnt} + ^ [x^/n\ - [ry/n\ 




n y/n 

P{B^ 2 t >x-r) = P{B 2 ^ <r-x) 


The first and the fourth sums in (5.26) are handled by Riemann sum arguments, with 
estimates to control the tails |i| > Cy/n, to obtain 

hm n-V2 y p^xlntnbiHrnV^} > 

2>0,z+/i;>0 

r*+oo 


>i + k') 


and 


n-\-oo 

/ ^ ^ - x)dx 

J 0 

Rb'i +Tji vn ^ ^ 


Urn n-':^ y p{xr';'"r"'‘' 

n—^■OD ' ^ ^ 


421 


2<0,2+fc<0 


< i + k) 


= / > rp - x)dx. 

J — OO 

More details can be found in [38]. By (5.16), this gives the conclusion. 


□ 
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We let re —)• oo in (5.20), with bounds (5.23) and (5.24) and with limit (5.25) apply the 
dominated convergence theorem, and recall the definition (3.1) of poj to conclude that 

(5.27) Ji^ ^1 = Qo (^2 > ^ 12 )) • 

i<ii J2 <a^ 

If the limit on the right of (5.27) vanishes then '^f=idjSn{tj,rj) converges weakly to 
zero. We continue by assuming that the limit on the right of (5.27) is strictly positive. 
The proof of Proposition 5.3 for case (b) will be completed by the following central limit 
theorem for linear processes due to Peligrad and Utev [26, Theorem 2.2(c)]: 


Theorem 5.6. Let {bn,i : 
(5.28) 

and 


—nin < i < run, re G be real numbers that satisfy 


lim sup b 

n—)-co 


2 

n,i 


< OO 


(5.29) 


lim max|6„ j| = 0, 
n^oo ’ 


where bn,i = 0 for |i| > m^. Let {z{i) : i £ 7,} be a centered, strongly mixing and non¬ 
degenerate (Var[z(0)] > 0) stationary sequence with strong mixing coefficients {a(j)} such 
that 


(5.30) 


Var 


m„ 

E 


bn,iZii) 


= 1 


and there exists 6 > 0 so that < 00 and + l)^^^Q:(j) < 00 . Then, 


(5.31) 


rrin 

bn,iz{i) ^ A7(0,1) as n ^ 00 . 

i=-ran 


Apply this to bn,i = an^ijon and z{i) = r]Q{i) — fj,Q. Since the random walk steps are 
bounded, an,i = 0 for large enough i when re is given. From (5.25) and (5.27), 



<00 


and condition (5.28) is satisfied. By (5.19) |a„_j| < re thereby 

1 ^ 

max|6„d < ——— > ^ 0 as re — 00 , 

iez ' ’ ' reV^dn 

verifying (5.29). Theorem 5.6 gives d~^ Siez *(%(*) “ To) =► A7(0,1) which, combined 
with (5.27), yields 

N 

y ^ QjSnifj ; ^j) ^ , go y ^ ((^ji J )) (^j2 5 ^12))^ ■ 

1 = 1 1<1iJ2<Y 

This concludes the proof of Proposition 5.3 of case (b). 
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Case (c). It remains to consider the case where the initial increment sequence 770 is 
dehned by (3.10). Let £(n) be any increasing sequence such that lim^^oo= 0. 
Split (5.18), with Ho = 0 now, into two terms: 

N £{n)-l 

= an,iC-fe(j)[p^(Lj) -/(* “ l>i)] 

(5 32) 

+ Y Y ^n,i^-kij)\p^ii,j) l,i)] = Tn,l+Tn,2- 

k=£{n) i,j& 

In the next lemma we show that the term is negligible. 


Lemma 5.7. With lim nj-Ji{n) = 0, we have lim E[|T„ 2 |^] 0. 

Proof. Due to the bounded range of the jump kernel p{0,j), after the Hrst equality sign 
below the sums over ii and 72 have 0(n) terms, and for a fixed k the sum over j is also 
finite. 

00 

E[ |r„, 2 p ] = X] Y (^n,han,i2 Y^PilJ “ Ph-l,j)ipi,j “ pYij) 

k=i{n) 11,12^"^ 

CX) 

~ Y1 Y1 ~ Qi2-ii+lfl ~ 1i2-ii-l,o') 

k=£{n) ii ,^2^^ 

00 

( 5 - 33 ) — Y1 Y1 ~ ^J+ 1.0 ~ ■ Y1 ^n,i(£n,i+j ■ 

jEZ k=i{n) i&t 

Bound the middle sum using the characteristic function (fq{9) = ■ 

00 

Yj “ ^j+ 1,0 ~ ^i- 1 , 0 ) 


k=i{n) 


1 r'^ 

S-r-L*-'"" 


e-ije _ ^-i{j+i)e _ ^-Hj-i)eyQ 


1 r 


1 - 4>q(e) 


<c (fY\e)de = Cq^^^\o,o) 


< 


c 




The ratio (1 — cos0)/(l — 4>q{9)) is bounded over [—vr, vr] because it has a finite limit at 
0 = 0 and by the span 1 property ( f ) q {0) = 1 only at 0 = 0. 

The sum over j on line (5.33) has 0{n) nonzero terms. Recalling (5.24), we conclude 

□ 


that E[ |T„, 2 p ] < Cn/ y^£(n) —^ 0 as n ^ 00 . 


Let Un,k{j) = i-k{j)Y.i&an,i{plj -Pt-i,j) so that T„,i = YYq ^T.j&Un,k{j) is a 
finite sum of independent mean zero random variables. To prove the goal (5.17), the 
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Lindeberg-Feller theorem shows that ^ '^f=i(^jSitj,rj). First the limiting variance. 

£(n)-l 

^ [2g^(j, 0) - q^(j + 1,0)- q^(j - 1, 0)] ^ an,ian,i+j 
j£l, k=0 iez 

(5.34) =^IY^ [«(•?' - 1) + “O' + 1) - 2a(j)] J]] an,ian,i+j - 

jez iez 

^ ^12^2 ( (^J1 ) ) {tj2 ) ^i2 )) • 

^ i<ii,i2<A^ 

The last limit came from a combination of (2.23) and Lemmas 5.5 and 5.7 . 

For the Lindeberg condition (condition (ii) of Theorem 3.4.5 on p. 129 in [10]) the task 
is to show that, for all e > 0, 


(5.35) ^ ^ E[f/2^(j)l{|C/n,fe(j)| > £}] ^ 0. 

jEZ k=0 

By Schwarz and Markov inequalities, 

i(n) — l i{n)—l 

E ^ui,{j)i{\Un,kU)\>^}] < ^E E ^[KM 

j£Z k=0 

(5.36) . 

jez k=o iez 

To two powers of the quantity in parentheses apply the inequality 

E [/(bj) +P^{i - 1, j)] < 

iez 


jSZ fc=0 

IE[Co"(0)] 


Then line (5.36) is bounded above by 

(5.37) ^ (E“->''0+')-/(^-l,j)]) 

jez k=o iez 

which vanishes as n —)• oo because the finite limit of the double sum of squares is exactly 
what was treated in (5.34). This verihes (5.35). Collecting the pieces, we have showed that 
7n,i This, together with (5.32) and Lemma 5.7, implies the goal (5.17). 

Proposition 5.3 has now been proved also for case (c) of Theorem 3.1. □ 


To complete the proof of Theorem 3.1, apply Propositions 5.2 and 5.3 to the right-hand 
side of decomposition (5.3), and use the independence of the processes Sn and 


6. Height fluctuations: process-level convergence 

In this section we prove Theorem 3.4. To simplify the exposition we take Q = [0,1]^. 
Theorem 2 in Bickel and Wichura [3] gives the following necessary and sufficient condition 
for weak convergence Xn X of i42-valued processes. 
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(i) (Convergence of finite-dimensional distributions.) For all finite sets {{ti,ri)}fLi C 
[0,1]^, we have 

{Xn{ti,ri),... ,Xn{tN,rN)) {X(ti,ri),.. .,X{tN,rN))- 

(ii) (Tightness.) Ve > 0, lim^^o P{t(;^(X„) > e} = 0, for the modulus 

t(;^(x) = inf max sup \x{t,r) — x{s,q)\, x € D 2 , 

^ GeA (t,r),(s, 5 )eG 

where the infimum is over partitions A of [0,1]^ formed by finitely many lines 
parallel to the coordinate axes and such that any element G of A is a left-closed, 
right-open rectangle with diameter at least 5. 

We proved the finite-dimensional marginal convergence (i) in Theorem 3.1. For the 
tightness proof that follows we check the sufficient conditions given by the next lemma 
(Proposition 2 in [21]) in terms of the modulus of continuity 

(6.1) ws{x) = sup \x{t,r) — x{s,q)\, x ^ D 2 . 

h,r),(s,q)e[ 0 ,l ]2 

|(t,r-)-(s,g )|<5 

Lemma 6.1. Let {Xn} he a sequence of D 2 -valued processes. Assume that there exists a 
decreasing sequence (5n \ 0 such that 

(i) 3 /3 > 0, A > 2, and C > 0 such that, for all large enough n, 

(6.2) E{\Xn{t,r)-Xn{s,q)f) < C\{t,r) - (s,g)r 
holds for all {t,r),{s,q) £ [0,1]^ at distance |(t,r) — {s,q)\ > 6n; 

(ii) Ve, r] > 0, there exists an uq > 0 such that for all n > no, 

(6.3) P{ws„{Xn) > e} < T]. 

Then, for all e,r] > 0, there exist 0 < 5 < 1 and integer uq < 00 such that 

P{ws{Xn) > e} < r] Vn > no. 

The two tightness conditions (i) and (ii) are checked in Lemmas 6.2 and 6.5. 

Lemma 6.2. Assume the assumptions of Theorem 3.4- Fix k and 7 such that 2 < k < 3 
and 0 < 7 < ^ and let 6n = n~'^. Then, there exists a constant C > 0 such that for all 
sufficiently large n, 

(6.4) E{\hnit,r) -hn{s,q)\^‘^) < C\{t,r) - (s,g)|'' 
for all t, s,r,q £ [0,1] with \ {t,r) — (s, g) j > n~'^. 

Proof. From decomposition (5.3), for a constant C, 

C-^E(}hn{t,r)-hnis,q)\^^) < \Hnit,r) -Hn{s,q)\^^ 

+E(i;F„(Lr) -F^s,^)]!^) +E(i:S„(Lr) -54s,g)li2) . 

Estimate (6.4) comes by treating each term on the right-hand side of (6.5) in turn. 
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As observed in the early part of Section 5, Hn{t,r) = 0{n ^1^) uniformly, and the first 
term on the right of (6.5) satisfies 


(6.6) \Hn{t,r) - Hn{s,q)\^‘^ <Cn ^ < Cn < C\{t,r) - {s,q)\'^. 


The following argument will be used more than once: if {Ci} are i.i.d. mean zero variables 
with a finite 12th moment and |ai| V 1 < il, then 


(6.7) 


E 

< C 


12 


^ ^ O-iCi j E 'y ^ ■ ■ ■ 0,ii2 Cii ' ' ' Ci_ 




E EKr Ei“.=rd- E 


I 


im/c 


l<fc<6 h 

mi>2: miH-|-f«fe=12 


12 


'^k 


where C is a combinatorial constant. 

For the second term on the right of (6.5), recalling (5.5), and with t> s, 


Fn{t,r) - Fn{s,q) 

lns\ 

= EE&w 


k=l xGT, 


+ n 


[nt\ 

k= [nsj +1 xGh 


= -) - = -) 
= x). 


Apply (6.7): 


E\Fn{t,r) - Fr,{s,qt^ 
[nsj 

< Cn-^ ' 

k=l x£Z 


1+ee ="■) - =") 

=i)^| 


[ntj 

E 

fc=[nsj+l icEZ 




■ [ntj —k 


6 
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The sums inside the braces develop as follows. 


( 6 . 8 ) 


[nsj 


k=l x£.'L 


\nth\ + \r^Jn\ 


\nt\-k ^\^[nsi-k ~ 


[ntj 

k=[ns\-^l x^Ta 


[nt\ 

p/YK6J + [rVriJ _ Y\ntb\ + \ry/n\ _ 
/ ^ \‘^\nt\—k [nt\—k 

k=l 


[ns\ 

p/j^[nsb] + [q^} _ j^\nsb\ + 'iqy/n\ _ 
/ V \ [ns\—k lns\—k 


k=l 


[nsj 


2 ^ + _ j^lnsb\ + lq^} _ „n 

/ V \ \_nt\—k \_ns\—k ’ 


k=l 


where X, and X. are independent random walks with transition probability p. 
Recalling transition probability q from (2.11), 


[nt\ 

'SX td( v-\ntb\ + \r^/n\ _ ^\ntb\ + \r^/n\ 
/ ^ ^ \^\r).t\—k \nt\—k 

k=l 


\nt\ [nt\ — l 

0 ) = J]gLntJ-fc(o^O)^ E 

k=l k=0 


Similarly, 


[nsj [nsj—1 

E = 0 ) = E 

k=l k=0 


And, 


[nsj [nsj—1 

p(^lntb\ + lry/n\ _ j^[nsb] + [q^\ _ p/j^lntb\ + lr^\ _ j^[nsbi + [q^\ 

/ ^ V [nt\—k [ns\—k ^ V [nt\ — [ns\-\-k k 

k=l k=0 


[nsj — 1 


E ® 




[ni6J + [rv^J 
[_nt J — L’T'SJ 


-X 


lnsb\ + lqy/n\ 


k=0 


0 ) 


= E 


[nsj — 1 

E fc / ^IntfeJ + Lr^/iil 


fc=0 


[nsh\ - [qV ^\, O) 
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In sum, we rewrite the right-hand side of (6.8) as 




[nsj — 1 


[nsj — 1 


^ 9^(0,0)+ q\0,0)-2El 


k=0 
[nt\ —1 


k=0 


k=0 


[nsJ —1 


(6.9) = ^ g"(0,0)+2^| ^ 

k= [nsJ 


k=0 


The first term in (6.9) is bounded by 

[nt\—l [nt\ — l 

( 6 . 10 ) 


[ntj-[nsj 

E E E —sc- 


Vk 


1 + y/ (t — s) 


n 


/c=[n5j fc=[72sj k = l 

To bound the second term in (6.9), observe first that the terms in the potential kernel (2.12) 
are nonnegative: 

q\x,Q) =Y,p\x,y)p\^,y) < \Y.[{p\x,y)f + {p\Q,y)f] =Y,{p\Q,y)f 

y& ysz 

= g"(o,o). 

Now the second term in (6.9) is bounded by 
( L’T’sJ-i 

E E [«*(».») -- i-'-j - L‘iVij,o)_ 


k=0 


<E 


(6.11) <C{\t — + \r - q\)y/n + C, 

where the second inequality is due to 

( 6 . 12 ) 


lim |x| ^a(a;) = 

from P28.4 on p. 345 of [34] and by the symmetry of a(x). 

Combining (6.10) and (6.11) gives a bound for (6.8), and hence the second term on the 
right of (6.5) satisfies 

(6.13) E (|F„(t,r) - F„(s,g)|^^) < Cn“^ + \r - q\^ y/n + 1 . 

For the third term on the right of (6.5), from (5.6), 

Snit,r) - Snis,q) 

= n-‘/->E(’'o(») - fo) [^(' s - p(i < 


i>0 


E(’)»(’) - w) [P(> > -P(i>^ 


. [ns6J + [qy/n\ > 
[nsJ > 


i<0 
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Note that = k + and define events 

^1,* = {^nti ^ - [ntb\ - < - [nsb\ - [gv^J} , 

A 2 ,i = {x[nti < - [ntb\ - > - [nsb\ - } • 

Then we can rewrite the difference above as 

(6.14) Sn{t, r) - Sn{s, q) = “ f^o) [P{^i,i) - P{^ 2 ,i)] ■ 

ieZ 


We prove an intermediate bound where the different cases (a), (b) and (c) of Theorem 3.4 
are felt. 


Lemma 6.3. Suppose the initial increment sequence {r/o(a:)} satisfies one of the assumptions 
(a), (b) and (c) of Theorem 3.f. Then we have the inequality 

(6.15) E[(54Lr)-54s,g))^"] <Cn-3{l + ^ 

Proof. Case (a). The case of i.i.d. {r]Q{x)} is handled by the argument in (6.7), beginning 
with (6.14). Details can be found in [38]. 

Case (b). For the case of strongly mixing initial increments we state a lemma from Rio’s 
lectures (see Theorem 2.2 and the derivation of equation (C.6) in [29]). 


Lemma 6.4. (a) Let m £ N, {Xjjigp} centered real random variables with strong mixing 
coefficients {a{k)}k>Q and define 

a~^{u) = inf{A: e Z+ : a{k) <u} = '^ !{«<«(*)}• 

i>0 


Assume < oo and let Qk{u) be the quantile function of \Xk\. Let Sn = ^k- 

Then there exist positive constants am and bm such that 


E {Sl^) < a, 


n 


(6.16) 


y^[a ^(u) A n]Ql{u)du \ 

V'^ k=i J 

n 

+ bm'^ [a~^{u) Anf'^~^Qf^{u)du. 
do 


(b) Suppose centered have finite rth moment. Then for p £ [l,r) there exists a 

constant Cp > 0 such that 

n „i 

/ [a~^{u) An]P~^Ql{u)du 

fc=i-^o 


< Cp ( ^(f + !)(■ 


1—p/r 


pr—2r-\-p) / (r—p) 


a{i)] j;(Eix,r 


\p/r 


A=0 


k=l 


(6.17) 
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Apply this lemma to representation (6.14). Note that G Z —P(A 2 ,i) 7 ^ 

0} = 0(n) due to the bounded support of p (2.1). Letting Qi be the quantile function of 
\r]o{-i) - Po\- I and m = 6 , (6.16) gives 


E[(54Lr)-5„(s,g))''] 


< Cn-^ 



a ^{u) /\kn]Qi{u)du 



Let p = 2,r = 


12 in (6.17), we can get an upper bound for 


/(I[a ^{u)Akn]Qiiu)du, 


VJ' [a ^{u) Akn]Qi{u)du 


iez 


5/6 


< C2 ( ) X] ' \P{Al,j) - P{^Aj)\ 


0=0 


jez 


/ \ 5/6 

< c Y.ii + [p(A,j) + p(A,j)] <cY;^ [p(Mj) + p(A^j)]. 

V*=o / jez jez 


By the same token, let p = 12, r = 12 + 5 in (6.17), we can show that 



/ fcn \ V(i2+5) 

< C £(i + l)i»+l32/«a(i) [p(Ai,,) + P{Aij)] 

\i=o / jez 


<CY.[P{Ai.l)+P(A^.l)\- 

iez 


Hence, 


(6.18) 


E[(5„(t,r) - 5n(s,g))^^] 

< Cn-4 [P(Ai,,) + P(A2,,)])' + [P(Ai,,) + P(A2 j)] 

^ jez jez 

< Cn 5|^1 + ^ ^ [H(Ai^m) + P(A 2 ,m)] j" • 

mSZ 


Case (c). From (6.14) and (3.10), 

00 

Sn{t,r) - Sn{s,q) = ^ - P0,j+i+l) [Pi^l,i) “ Pi^2,i)] ■ 

i,jG'^ k=0 
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Then by (6.7) 


E[(5n(t,r) - Sn{s,q)) 


12n 


(6.19) 


2 \ 6 


< Cn-'‘ ( 1 + E E E (rf. j+. - pj. J+<+i) [P(-4 m) - P(-4w)] 

^ jGl^k=0 

Expand the sum of squares inside the parentheses: 

oo 

EE E ^Po,j+ii Po, j+ii+l)(7*0, j+i2 Po, j+i^ + l) 

k=0 ^1,22^=^ 

X [P{A,,i,) - P{A2,i,)] [P{Ar,i,) - P{A2,,)] 

= W - 1 ) + «(^ + 1 ) - M^)] • [^(^ 1.0 - ^(^ 2 , 0 ] [P{^i,i+e) - P{A2,i+e)] 


< [a{i — 1) + a{l + 1) — 2a{l)] 

£ez 

^ T.^{lP(Av)-P(A2.P+[P(Ai,t+i)-P(A,,+,p} 

iez 

=iE - p(A,)]' < 4 e +-p(-42,i)]. 

iez iez 

where the last equality is from (2.23). Inserting this into (6.19) gives 

E [(5„(t, r) -Sn{s, q)) < Cn-3{ 1 + [P(Ai^,) + ^(Tla,*)] 

iez 


The proof of Lemma 6.3 is complete. 

We continue by bounding the sum on the right of (6.15). Suppose t > s. 

P{Ai^m) = Y Pi^fntl ^ - [ntb\ - [rVnJ - < - [nsb\ - [q^/n\ - m) 

va£.'L mEZ 

= E E 

mSZ £>m 

X > [nsJiJ - lntb\ + [q^\ - [r^/ij -m + tj 

= — [ns6J — yq\/n\ — k — m) 

ttiEZ fc>0 

X P{Ynt\-lns\ > lnsb\ - [ntb\ + [gv^J - lrVn\ + k) 

(6.20) =YP{Ynti-[nsi > lnsb\ - [ntb] + lqy/n\ - [ry^J + k) 

k>0 

Similarly, 

(6.21) Y P{^2,m) =Y^i^lnt\-lns\ < + ^) ' 
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Combining (6.20) and (6.21), 


( 6 . 22 ) 


= '^^i^lnt\-lns\ - lnsb\ + [ntb\ - [q^\ + [rV^J > k) 
k>0 

+ - l'^sb\ + [ntb\ - [qV^\ + [rv^J < k) 

k<0 


= '^P{\^[nt\-[ns\ - lnsb\ + [ntb\ - [q^/n\ + | > k) 

k>0 

< - [ns 6 J + [ntb\ I + |r - g|v^ + 1 

< C[^/{t - s)n + |r - q\\^ + l). 

Combining (6.15) and (6.22) gives this bound for the third term in (6.5): 

(6.23) E (|S'„(t,r) - 5n(s,g)|^^) < Cn~^[{y/\t - s\ + |r - q\)Vn + l]®. 

Return to (6.5) and apply (6.13), (6.23) and ( 6 . 6 ) to conclude that 

E[|7i„(t,r) -hnis,q)\^‘^] < Cn~^[{y/\t - s\ + |r - q\)Vn + l]® 

<C{\t-s\^ + \r-qf + n-^) < C (|t - sr + |r - g|''). 

In the last step we used t, s, r, g G [0,1], 2 < k < 3, and n~^ < = 6'^ < |(t, r) — (s, q)\^. 

This completes the proof of Lemma 6.2. □ 

The second tightness condition is verified as follows. 

Lemma 6.5. Under the conditions of Lemma 6.2, for any fixed 1 < 7 < 3/2, for Me > 0, 


lim P<^ sup |/in(t, r) -/in(s,g)l > e > = 0. 

I (7r),(s,q)e[0,l]= J 

|(7r)-(s,Q)|<n -< 

Proof. Dehne the intervals I{k) = [{k — l)n“^, (k + l)n“'>'] C [0,1]. For hxed e > 0, hrst by 
a union bound and then by decomposition (5.3), 


(6.24) 

(6.25) 

(6.26) 


sup \hn{t,r) - hnis,q)\ > £ 

0,r),(s,q)e[0,l]2 ^ 

\{t,r)-{s,q)\<n—y 

[n~<\ [n^J 


< V'P| sup \yLQHn{t,r) - iXQHn{kin ,k 2 n '^)| > ^ | 
tei{ki),rei{k2) 0 J 


ki=l k2=l 

[n<\ [nT-J 

+ 


y] y] P| sup \Sn{t,r) - Sn{kin ,k 2 n '^)| > || 

k,=lk,=l ^ tmki),r&i{k2) hJ 

[n"'J [nT-J 

+ y] y] p| sup \Fnit,r) - Fn{kin~'^ ,k 2 n~'^)\>^Y 


k\=l ^2=1 


te/(fci),re/(fc2) 
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Line (6.24) vanishes for large n because Hn{t,r) = 

For the second and third sums (6.25)-(6.26), observe from (5.6) and (5.5) that Sn{t,r) 
and Fn{t,r) depend on their argument {t,r) only through [nt\, [ntb\ and [r-y/nj. For large 
enough n, and any t G r G I{k 2 ), 

\nt — kin^~'^\ = n\t — kin~'^\ < n^~'^ < 1 / 2 , 

\nth — kiTn}~'^h\ = n\b\ ■ \t — kin~'^\ < \b\n^~'^ < 1/2, 

\r^/n — k2n^^‘^~'^\ = — k2n~'^\ < < 1 / 2 . 

Thus for t G I{ki) and r G I{k 2 ), each of [ntj, [nt 6 J and [r^/n\ can have at most one 
jump. For example, [ntj can only jump from [kin^~^\ — 1 to or from 

to [kin^~'^\ +1. As a result, 5„(t,r) and Fn{t,r) can take only at most 8 different values 
on I{ki) X I{k 2 )- 

Suppose {^^(L, rj)}i<j<£ with i < 8 and {ti,ri) G I{ki) x 1 (^ 2 ) captures the different 
values of Snit,r) on the square I{ki) x I{k 2 )- Then, 

p| sup \Sn{t,r) - Snikin~'^,k 2 n~'^)\>^ 
lte/(fci),re/(fc 2 ) ® 

e 

< ^P{|S'„(L,ri) - Snikin~'^,k 2 n~'^)\ > e/ 6 } 

e _ g 

< Cn~^e~^'^ ^ - kin~'^\ + {n - k2n~'^\J ^/n + 1 

2=1 

< + n“^)y/n + l]® < Cn~^ 

where the second inequality comes from Markov’s inequality and (6.23), and the third one 
from 7 > 1 . 

Therefore, for any 1 < 7 < 3/2, as n —>■ 00 , 

[nTJ \n<\ f \ 

Y Y p 1 sup |5„(i,r) - Sn{kin~^,k 2 n~^^)\ > | 1 < ^ 0. 

rel{k2) 

Same reasoning with inequality (6.13) shows that line (6.26) vanishes as n —>■ oo. Lemma 
6.5 has been proved. □ 

The proof of Theorem 3.4 is complete. 

Appendix A. An ergodic lemma 

In this section p{x, y) = p(0, y — x) is an arbitrary random walk kernel on with the 
property that the smallest subgroup that contains the support {x : p( 0 , x) > 0 } is itself. 

Lemma A.l. Transition probability p{x,y) has no bounded harmonic functions other than 
constants. 
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Proof. Suppose /i is a bounded harmonic function on that is, h{x) = YlyP{^>y)h{y) 
for all X £ Suppose p{u, v) > 0. Then the coupling described on p. 69 of [23] works 
to show that h{u) = h{v). (This coupling is also spelled out in Section 1.5 of [31].) By 
assumption for any G Z"^ there exists a path x = xq,xi, ... ,Xm = y in Z*^ such that 
p{xi, Xj+i) + p{xi+i,Xi) > 0 for each i = 0,... ,m — 1, and thereby h{x) = h{xo) = • • • = 
h{xm) = h{y). □ 

Let 5 be a Polish space, T = 5^“^, and shifts {0xC){y) = Ci^ + y) for G T and x,y G Z'^. 

Lemma A.2. Let v be a probability measure on T that is invariant and ergodic under the 
shift group and f G L^(p) with finite mean E'^[f]. For x G Z*^, t G Z+, and C £ L, define 

(A.l) gt{xX) = '^p\x,y)f{9yC). 

yeZ 


Then Vx G Z'^, gt{x, C) —5- E'^[f] in L^(p) as t —)■ oo. 

Proof. The argument is a Fourier analytic one suggested by the proof on p. 30-31 of [24]. 
The characteristic function of the jump probability is 

4>x{a) = '^p{0,y)e‘‘‘^'y, a£R^. 

yeZ 


For 0 < r < oo define truncated functions /r(C) = (“^) V (/(C) A r) and gt{x,r,() = 
X^ygzP*(a;,y)/r( 6 *yC)- Note that gt{x,r,C) is bounded, uniformly over x and C- We show 
that, as t ^ oo, gt{x,r,C) converges in L‘^{v) to a constant. An approximation via the 
truncation then implies the resnlt. 

The function ^(''^(x) = E'^[fr{C,)fr{0xC,)] is nonnegative definite (i.e. Ylixy^^''\^ ~ 
y)zx'^ > 0 for any choice of finitely many complex numbers {zx}). By Herglotz’ Theo¬ 
rem (Chapter XIX .6 in Feller [11]), there exists a bounded measure 7 on [— 7 r, 7 r)‘^ such 
that 

Let Xt and Xt be i.i.d. copies of the random walk with transition p. Compute: 

j 9t{x,r,C)gsix,r,C)i^{dC) = J E^[fr{6xtOW[fr{0x^C)]i^{dC) 

= I MexX)fri0x,CHdC) = E^^’^^[v^^\Xs-Xt)] 

= J [4>x{a)Y[f’x{a)Y'y{da) = J [4>x{a)Y[4>xia)Y'y{da). 
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The last equality came by switching s and t around in the calculation. From this, 


j [9t{x,rX) - gsix,rX)]^i^{dC) 

= j [9tix,r,C? - 2gtix,rX)gsix,rX) + 9six,rX)‘^]i'{dC) 

t>x{a)\‘^^ - 24>x{ay4>x{aY + \(l)x{a)f^]'y{da) = j \(j)x{aY - (j)x{ay\'^ l{da) 


(A.2) 


[ \(px{aY - (l>x{ayf lida). 

Ja^O 


Since is the smallest subgroup that contains the support of p, |(/>j\:(a)| < 1 Va G [—vr, 7r)'^\ 
{0} (T7.1 in [34]). Thus the integrand in (A.2) is bounded and converges to zero as 
oo. Hence {gt{x,rX)}te‘Z+ is Cauchy in Lp‘{v) and 3 a bounded L‘^{v) limit g(x,rX) = 
limt^^gt{x,r,C)- 

Letting s oo in gs+t{x,r,C) = Y.yP^{x,y)gs{y,r,0 implies 


g{x,r,C) = '^p\x,y)g{y,rX)- 

y 


Hence for almost every C, g{- ,rX) is a bounded harmonic function for p{x,y) and thereby 
a constant in x. Combining this with a shift gives g{0,r,() = g{x,r,() = 5(0, r, 0^:0 
Vx G By ergodicity, g{0,rX) equals i/-almost surely a constant, and hence g{x,r,C) 
equals i^-almost surely the same constant for all x G Z*^. 

Now we transfer these properties to gt{x, C) through approximation. By shift-invariance 


\\gt{x,C) - 9tix,r,C)\\i <'^pyx,y)\\f oOy - froOyWi = \\f - frWi —^ 0. 

y 

From this we deduce that {gt{x, •)} is a Cauchy sequence in L^X). Hence we have an L^X) 
limit g{x, C) = lim^^oo gt{x, C) and we can take t —)• oo in the bound above to get 


ll5(a:,C) -9{x,rX)\\i < ||/ - /r||i- 


Letting r —)• oo takes the right-hand side to zero, and we conclude that g{x, Q is p-almost 
surely a constant independent of x. This constant must equal E'^lf] because by the L^X) 
convergence E'^[g{x, •)] = limt^oo E’'[gt{x, •)] = E‘'[f]. □ 
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